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We consider cosmological models in which dark matter feels a fifth force mediated by the dark energy scalar 
field, also known as coupled dark energy. Our interest resides in estimating forecasts for future surveys like 
Euclid when we take into account non-linear effects, relying on new fitting functions that reproduce the non¬ 
linear matter power spectrum obtained from N-body simulations. 

We obtain fitting functions for models in which the dark matter-dark energy coupling is constant. Their 
validity is demonstrated for all available simulations in the redshift range z = 0 — 1.6 and wave modes below 
k = lOh/Mpc. These fitting formulas can be used to test the predictions of the model in the non-linear regime 
without the need for additional computing-intensive N-body simulations. We then use these fitting functions 
to perform forecasts on the constraining power that future galaxy-redshift surveys like Euclid will have on the 
coupling parameter, using the Eisher matrix method for galaxy clustering (GC) and weak lensing (WL). We find 
that by using information in the non-linear power spectrum, and combining the GC and WL probes, we can 
constrain the dark matter-dark energy coupling constant squared, j3^, with precision smaller than 4% and all 
other cosmological parameters better than 1%, which is a considerable improvement of more than an order of 
magnitude compared to corresponding linear power spectrum forecasts with the same survey specifications. 


I. INTRODUCTION 


Cosmological data have reached a level of precision that forces us to address effects which may have been safely neglected 
in the past. Among them, predicting the behaviour of theoretical models at non-linear scales gives a more concrete chance 
to disentangle ACDM from alternative scenarios. The non-linear power spectrum contains useful information encoded in the 
clustering of galaxies and the weak lensing correlations, which can help to constrain more tightly the cosmological parameters 
and remove certain degeneracies between them. The drawback relies in the fact that the full estimate of the non-linear power 
spectrum requires time-demanding and computationally expensive N-body simulations. As a consequence, there exist simu¬ 
lations only for a handful of models beyond ACDM on a very limited subset of the parameter space of dark energy theories. 
This makes it hard to use the non-linear power spectrum in forecasts (see e.g. CD, both when doing cosmological parameter 
estimation via Monte Carlo simulations and with the simpler approximation of the Fisher matrix formalism. The alternative is 
to employ fitting functions for the non-linear power spectrum, like Halofit Hal ID or machine learning estimators like the cosmo¬ 
logical emulator (CosmicEmu) ||4]-0 and the interpolator PkANN 17] [H based on artificial neural networks. All these methods 
are calibrated or trained using large sets of N-body simulations and typically permit an exploration of the parameter space just 
around a fiducial ACDM cosmology or wCDM cosmology with a constant w. An extension of Halofit to account for massive 
neutrinos was realised in Q, which was then used in IfTOl to perform a parameter estimation using a Markov Chain Monte Carlo 
technique. Other authors have recently used polynomial fits to take into account systematic effects in the non-linear regime like 
baryonic effects such as in im or to parametrize the cosmological dependence of non-linear clustering, beyond the Zeldovich 
approximation (see El). These fitting functions have, however, an intrinsic error that usually increases when scales become 
highly non-linear; therefore one has to be aware of the range in scales and redshifts they were designed to work on, in order to 
keep the errors induced on the power spectrum under control. 

As an alternative to fitting formulae and N-body simulations, there has been in recent years a substantial progress in semi- 
analytical methods to calculate the non-linear power spectrum at mildly non-linear scales, such as for example renormalized 
perturbation theory El and all other resummation methods derived from it (for example lfT4ll L the time renormalisation group 
(TRG) IITSi and effective field theories of large scale structure IfT^ . Most of these approaches claim to reach percent accuracy 
at the baryon acoustic oscillation (BAO) scale and a bit beyond the second BAO peak, but they are still not able to predict what 
happens at highly non-linear scales, when the single stream approximation does not hold any longer. 

In this paper, we will use N-Body simulations to find fitting functions for a class of models beyond ACDM usually refered 
to as coupled Dark Energy (CDE). These models, widely discussed in literature ifTTl - EOl . involve an extra degree of freedom, 
associated to a scalar field that provides acceleration and mediates a fifth force, in addition to gravity, which is felt by dark matter 
particles only. Semi-analytical non-linear analysis lITTHTTll and cosmological N-body simulations within coupled Dark Energy 
have been performed by many different groups Il23l428]l and their effects on large scale structure formation have been identified 
and characterised. The power spectrum, halo mass functions and concentration, halo spin and sphericity, voids and amount of 
substructures show noticeable differences compared to a simple ACDM model (see for example Il29l - l3^ and the review article 
03). Eor a constant coupling, constraints have been found for a variety of probes ll36]439l . with the latest ones discussed by the 
Planck collaboration in BOl . Recently there have been attempts to constrain more general couplings between dark matter (DM) 
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and dark energy (DE) using large-scale structure ED, CMB l42l or laboratory experiments ll43l . Forecasts on coupled dark 
energy using galaxy clustering (GC) and weak lensing (WL) measurements for future surveys like Euclid have been discussed 
in P4ll . but have been performed using only linear power spectra for the CDE models. The TRG method has been extended to 
coupled Dark Energy Il45]l and to massive neutrinos ll46l . However the TRG method does not produce a reliable estimation of 
the power spectrum for scales larger than k « 0.3/i/Mpc, which makes them less suitable for forecasts which attempt to extract 
information on highly non-linear scales. 

The Codecs (Coupled Dark Energy Cosmological Simulations) set of N-body simulations Il27ll has shown that CDE models 
have characteristic and measurable features in the morphology and history of non-linear structures, such as halos, subhalos and 
voids, and therefore in the non-linear power spectrum. 

The aim of this paper is to create fitting functions which are valid in the observable regime of non-linear perturbations at all 
interesting redshifts and reproduce the subtle effects of coupled dark energy on the non-linear power spectrum while allowing 
us to vary the different parameters of the model. We use this to perform forecasts of cosmological parameters assuming coupled 
dark energy as the fiducial model, using galaxy clustering and weak lensing as observational tools, as expected for future surveys 
like Euclid Il47ll48l . We do a careful treatment of errors and systematics, so that we take into account all errors induced by our 
fitting functions, the cosmic emulators and the extraction of the power spectrum from the N-body simulation into the analysis, 
systematics related to non-linear effects that affect the redshift space distortions and the lensing signals. In this way we obtain a 
conservative estimate on how well a probe like Euclid, will be able to measure a DM-DE coupling. 


II. COUPLED DARK ENERGY 

In this section we briefly review the main equations governing the coupled DE models that will be investigated in the present 
work. For a more detailed and extended discussion and for the full derivation of the equations we refer the interested reader to 
some of the original publications on the subject lfT8l - l20l 1^ . 

The background evolution for the coupled DE scenario model is described by the following equations, in which the subscripts 
r, b, c and (j), indicate radiation, baryons, cold dark matter (CDM) and the dark energy scalar field, respectively: 


Pc + 3//p. = -y|/3(0)^, (2) 

Pb + 3Hpi, = 0, (3) 

Pr+4Hpr = 0 , (4) 

~ ~rjTiPb3-Pc +Pr +P(j) ) ■ (5) 

Mpi 


We express from now on the scalar field p in units of the Planck mass Mpi = 1 and choose as potential V{p) an 

exponential V{p) = Il49ll50l . The coupling function P{p) defines the strength of the interaction between the DE fluid 

and CDM particles and in the present work we will restrict our analysis to the simplified case of a constant coupling (0) = j3, 
although in general it could be a field-dependent quantity IHIIED. 

The current constraints on a coupling to ordinary matter are very tight. The “post-Einstein” coupling parameter 7 that measures 
the local admixture of a spin-0 field to gravity is constrained in Solar System experiments (see e.g. the PDG review 15^ and also 
ED El) roughly to |7| < 4 • 10 ^ . This parameter enters the modification of the effective Newton constant as Ge/f — Gn{1 — 
7/2); in our notation, this is Gg// = GAf(l +4p^/3) (see below) and therefore = —37/8. A coupling Plaryons appears then 
constrained to be smaller than 10^^ roughly, and we assume therefore that is completely negligible. As a consequence, baryons 
follow the usual geodesics of a FLRW cosmology, which allows coupled DE to pass the stringent local gravity constraints 
without the need to employ any screening mechanism EH- 

Due to the exchange of energy between DE and CDM, the energy density of the latter will no longer scale as the cosmic 
volume, and by assuming the conservation of the CDM particle number one can derive the time evolution of the CDM particle 
mass by integrating Eq.Q between the present time (z = 0 ) and any other redshift z: 

mc{z) = . ( 6 ) 

At the level of linear perturbations, coupled DE models are characterised by a different evolution of the baryonic and CDM 
density fluctuations, as a consequence of the selective interaction between DE and CDM particles only. In the sub-horizon limit. 
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Parameter 

Explanation 

Value 

Reference ACDM 

A 

Potential normalization 

0.0218 

- 

a 

Potential slope 

0.08 

- 

II 

Scalar field normalization 

0 

- 

p 

Coupling parameter 

{0.05,0.10,0.15} 

- 

(z = 0) 

DE equation of state { 

-0.997, -0.995, -0.992} 

-1 

CT8(z = 0) 

Power spectrum amplitude 

{0.825, 0.875, 0.967} 

0.809 


Table I. The main parameters and normalizations of the three coupled DE models of the CODECS suite considered in the present work. The 
last column displays the corresponding values for the reference ACDM cosmology. 


for which aH/k <C 1, linear perturbations in coupled dark energy follow the equations lfT9ll20ll : 


L = -2H 


l-j3 


0 


hV6 

5h = —2H 5h + AkG [Pb^b + Pc^i 


5c + AkG [Pb5b + Pc^cF c 


(7) 

( 8 ) 


where = 1 + represents the effective “fifth force” acting on the CDM particles. The term proportional to in equation 
0 is a velocity-dependent term that modifies the standard cosmological friction; this arises as a consequence of momentum 
conservation for the CDM particles and has a considerable effect on structure formation ll2^ 12411311 . Since baryons are uncou¬ 
pled, their perturbations evolve according to the standard equation. Nonetheless, baryons will still be indirectly affected by the 
coupling as the source term on the right-hand side of equation (|^ includes the CDM density perturbations. 

At the level of non-linear perturbations, several methods have been devised to predict the small scale effects of coupled 
DE, from semi-analytical methods like spherical collapse ll5^ . to time renormalisation group ll22l to full N-body simulations 
ll^ 1^ [57] [5^ . In the following, we will work with the publicly available data of the CODECS simulations (27) that represents the 
largest set of cosmological N-body simulations for coupled DE models to date. 


III. THE CODECS SIMULATIONS 


The Codecs suit^lncludes simulated periodic volumes of the universe at different scales and with different physical ingredients (as e.g. 
simulations with and without hydrodynamics) in the context of a series of coupled DE cosmologies characterised by various choices of the 
self-interaction potential V{(j>) and of the coupling functions jS(p). The simulations have been performed with a suitably modified version 
(see {2l\ for more details on the numerical implementation) of the widely-used TreePM N-body code GADGET Il5^ . Such modified version 
includes all the relevant effects that characterise coupled DE cosmologies, from the modified background evolution to the CDM particle mass 
variation, the “fifth-force" and the velocity-dependent acceleration appearing in equation 0. 

As already stated above, in this work we will consider - besides the reference ACDM simulation - the subset of CODECS 
runs characterised by an exponential potential and by a constant coupling function. This consists of three different coupled DE models with 
the same potential slope a and with three values of the coupling j6 = {0.05,0,1,0.15}. The short names for these simulations are respectively 
EXPOOl, EXP002 and EXP003. All the models have been built in order to have the same cosmological parameters at z = 0 consistent with the 
WMAP7 results (60), see Table jn] with the obvious exception of the value of the equation of state parameter wq, that changes from model to 
model due to the different dynamics of the DE scalar field. The present observational constraints on the cosmological parameters have only 
slightly changed with the latest updated release of Planck data ED with respect to the assumed WMAP7 values, and are still good enough 
for the purposes of this paper, being in good agreement with large scale structure observations. Eor what concerns linear perturbations, all 
cosmologies have been normalised to have the same statistics (i.e. the same power spectrum shape and amplitude) of density fluctuations at the 
redshift of the Cosmic Microwave Background zcmb ~ 1100. As a consequence of this choice and of the different growth associated with the 
various coupling values, all the models will have a different normalisation cfg of the linear perturbations amplitude at z = 0. The main features 
of these models are summarised in Table jn] We refer to (23 for further details. 

For the purposes of the present work, we will employ the matter power spectra extracted from the CoDECS runs of these 
cosmologies at different redshifts in order to find a fitting formula that captures with high accuracy the deviations of the coupled 
DE nonlinear power spectra from the reference ACDM case. 


* see also the public CoDECS database at www.marcobaldi.it/CoDECS/ 
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Parameter 

Value 

Ho 

70.3kms^*Mpc^* 

^CDM 

0.226 

^DE 

0.729 


2.42 X 10^® 

Slh 

0.0451 

ns 

0.966 


Table II. The set of cosmological parameters used in all CoDECS simulations, consistent with the WMAP7 results. 


IV. EXTRACTING THE POWER SPECTRUM AT SMALL SCALES 

The power spectrum P{k), is defined as the ensemble average of the density contrast in Fourier space {5{k)5{k')) = 
{27z)^P{k)dD{k + k'). When trying to extract this quantity from an N-body simulation, one has to take into account several 
technicalities related to sampling effects which appear as a consequence of treating a discrete distribution of particles. For more 
details related to this problem and different solution methods, see ||62l. 

At large scales (or equivalently small k) the power spectrum suffers from uncertainties due to the finite size of the simulation 
box, since there are only few independent modes to sample the signal from. On the other hand, at high k, one is limited by the 
resolution of the simulation, since one cannot sample wave modes smaller than the typical grid size L/N where L is the length 
of one side of the simulation box and N is the number of particles. This maximum frequency is called the Nyquist frequency 
kf^y = 2kN/L and modes smaller than this cannot be reliably measured (this corresponds to the so-called aliasing effect). For 
the Codecs simulations used in this work, the Nyquist frequency has the value kjvy « 2.2 h/Mpc at present time. 

The power spectrum computation embedded in GADGET-3 that was adopted in the CODECS simulations employs the so-called 
folding method developed by l6^ - which is based on (Ml - to calculate the matter power spectrum for smaller scales than the Nyquist 
frequency. Following this method one ends up with two separate power spectra, P/op which is calculated using the simulation particle mesh 
(PM) at k < kj\fy and Pfghj which is the folded power valid for k > In order to provide a single sampling of the power spectrum across kffy 
the Codecs project employed a simple interpolation procedure around k^y by averaging the two power spectra in the region of 
overlap. Flowever, this might introduce some spurious features that appear only when considering the ratio of two power spectra 
P{k) at highly non-linear scales. Although such features are harmless for most practical purposes, for the aims of the present 
work it is very important to have accurate ratios of power spectra, since we want to calculate fitting functions that can capture the 
effect of the dark energy - dark matter coupling compared to ACDM and therefore we need to correct for these spurious effects. 

We then developed a PYTHON code that finds the optimal interpolation and matching between the Pf^p and Pfoid power 
spectra. By evaluating the first and second derivatives of the ratios and minimizing abrupt changes, it finds the optimal number 
of points to be removed from Prop due to the aliasing error and the number of points to be removed from the Pfoid due to the low 
sampling error; at the same time it looks for the optimal linear interpolation weights between them. Moreover, it cuts off the 
power spectrum when the shot noise error (Pshot = 1 /N) reaches 10% of the estimated power spectrum. This method improves 
considerably the convergence of the fitting functions at non-linear scales, allowing us to reach our accuracy goal of 1-2% (see 
section|v]on fitting functions). In figure Q, the uncertainties on P{k) are plotted, and the shaded region represents the error on 
the power spectrum due to the finite number of modes available. The clear jump present in this shaded region occurs at the scale 
in which the folded and top level power spectra have been matched together, which corresponds roughly to kj^y. 


V. FITTING FUNCTIONS 

A. The net effect of the DM-DE coupling at non-linear scales 

We now model the net effect of a coupled DE model with an exponential potential and a constant coupling on top of a fiducial 
ACDM non-linear power spectrum, by evaluating the ratio of the non-linear power spectrum of the coupled DE model, with 
respect to the one in the ACDM model, both extracted from the CoDECS simulations. In figure Q we show the ratio R{k\ j3 ,z) = 
PExp{k',P,z)/PLCDM{k',l5,z) as a function of the scale k, for four different redshifts z = {0, 0.55, 1.0, 1.61}; each panel contains 
the curves corresponding to the three available constant couplings in the CoDECS simulations: j3 = {0.05, 0.10, 0.15}. Since all 
coupled DE models have the same amplitude of perturbations at recombination, an increasing coupling has the effect of inducing 
a higher linear normalization ffg of the power spectrum at z = 0 m- Therefore, in order to see the net effect of the coupling 
at non-linear scales, each power spectrum ratio has been re-normalized by dividing each model by its respective so that 
at linear scales k < O.lh/Mpc all the ratios are unity. The net effect of the fifth force is a “bump” at non-linear scales, whose 
amplitude increases with higher couplings and whose maximum is shifted into higher wavenumbers k for higher redshifts. This 
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k [h/Mpc] k [h/Mpc] 


Figure 1. Ratios of the non-linear matter power spectra of the CODECS CDE models with an exponential potential, with respect to the 
Codecs ACDM power spectra, normalized by their respective CTg, evaluated at four different redshifts and three different coupling constants 
P. Upper left panel: z = 0, upper right panel: z = 0.55, lower left panel: z = 1.0, lower right panel: z = 1.61. The blue line represents the 
model with strongest coupling j3 = 0.15, the yellow line j3 = 0.10 and the red line the smallest available coupling j3 = 0.05. In order to be able 
to observe the net effect of the j8 coupling at non-linear scales we have divided each power spectrum by its respective CTg . 


extra information imprinted into the non-linear power spectrum is what we want to use to improve the estimation of parameters 
using future surveys. To achieve this, we will fit these curves which are functions of redshift, physical scale and coupling, using 
the minimal number of numerical parameters possible, while keeping the accuracy goal at the 1% level. 


B. Generating the fitting functions 

The fact that one can observe a clear trend that relates the amplitude of the signal to an increase of the coupling, together with 
a shift of the peak towards larger length scales when looking at smaller redshifts, is an indication that we should be able to hnd 
a relatively simple fitting formula describing this behaviour, which will be then a function of z, k and the coupling constant j3 
only. 

To perform the ht, we use a least-squares-minimizing technique, using the conjugate gradient method ll65ll . Taking into 
account the particular form in k-space of the ratios /?(k;j3,z) = PEXp{k',P,z)/PLCDM{k',P,z) that we need to ht, we use as an 
Ansatz different sigmoid functions to reproduce the particular form of the peak. For each htting model, we keep the same 
functional form for the k-dependence of /?(k;j3,z) at all redshifts and for all couplings. We tried 7 different sigmoid functions 
as htting models, but we only show the best two models M2 and M7 in table m All models contain 5 coefficients, which are 
dependent on the coupling j3 and the redshift z: (a,- = ao, ai, c, b, ko, with i — 0, ...,5). The coefficients ko and b determine 
qualitatively the form of the peak to be htted, while the others control mostly the shifting and the battening of the peak. 

Each coefficient a,- is then htted using a polynomial in j3 and z, up to a maximum of third order in powers of p and z. 
Polynomials of order 4 and 5 were also examined, but the gain in goodness of ht was minimal compared to the increase in the 
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Model Name 

Functional form 

# of coefficients 

R^-value 

M2 

f{k) = 1 +ao + ai ■ k + c ■ k ■ ~ Icq) ■ b) 

5 

0.99996 

M7 


5 

0.999989 


Table III. Fitting models M2 and M7 with their corresponding number of fitting coefficients and their -value. Each coefficient ao, ai, c, b 
and fco is fitted as a polynomial in the coupling parameter j8 and the redshift z (see Appendix|^for further details). The R^-value is a measure 
of the goodness of fit: a value of 1 corresponds to a perfect fit, while 0 means that the model does not fit the data. 


number of free parameters. Therefore, third order polynomials were the best compromise between complexity and goodness of 
fit. 


The best fitting models were chosen according to their coefficient of determination (also known as R^-value), which is a 
statistical measure for the goodness-of-fit Il66ll . It can be simply defined as = 1 —Sres/S,ot, where Sres is the residual sum 
of squares (the residual between the data points and the fitting function) and Stot, which is the total sum of squares and is 
proportional to the variance of the data. An R^-value of 1 corresponds to a perfect fit. The analytical expressions for the best 
models M2 and M7 are shown in table III together with their R^-value. We performed the whole Fisher analysis (see section 
|VII| below), for both of the models and the results on the parameter estimation are basically the same (less than half of a percent 
relative difference in the estimated final errors). 


C. Fitting functions and cosmological parameters 

In order to use the fitting formulae obtained before to forecast cosmological parameters of the model, we need to assume that 
the shape of the non-linear coupled DE signal does not change dramatically if the other cosmological parameters, apart from j3 
and (Jg, are modified by small amounts. This is justified since in the deeply non-linear regime, the evolution of perturbations 
is ruled by mode-mode coupling between high k wavevectors (non-linear k modes are not independent of each other anymore), 
which erases most of the information about initial conditions and makes the shape of the non-linear power spectrum at large 
k practically independent of cosmological background parameters, such as tig, and h. This has been shown to be the case 
when calculating perturbatively non-linear corrections to the power spectrum, see for example ||67]|68], but also analyzing the 
covariance matrix of non-linear power spectra using a large suite of N-body simulations as was investigated in ll69l . Furthermore, 
the decoupling of virialized structures in the small scale regime from the background dynamics of the universe, is one of the 
cornerstones of the recently developed effective field theory of large-scale structure ifTOl and was also shown to be approximately 
true using a coarse grained cosmological perturbation approach ini. Since j3 and ffg are quantities that affect directly the linear 
perturbations and the virialization dynamics and we are only looking at a particular signal at very small scales, they should be 
the main parameters determining the shape of the non-linear “bump”. An investigation of how robust these fitting formulas are, 
with respect to a change of cosmological parameters, would need either more high-resolution N-body simulation for coupled 
DE scenarios or the development of a consistent perturbation theory for modified gravity and scalar-tensor theories that reaches 
highly non-linear scales. 


VI. NON-LINEAR POWER SPECTRUM AND ERROR ESTIMATION 

The accuracy of the fitting functions when compared to the original N-body simulation power spectra is shown in figure]^ 
where the absolute value of the percentage error between fitting function and the original power spectra is plotted as a function 
of redshift and scale. In this case we show that for the model M7 the etTor remains well below 1% for a coupling constant of 
P — 0.05 and for the scales and redshifts we are interested in. Eor higher couplings the etTor goes up to a maximum of 1.5%. 
Eor the model M2 the same test gives similar results, yielding <1% etTor in all interesting scales and redshifts. We can then 
expect that when applying these fitting functions on top of ACDM power spectra, the intrinsic etTor is less than 1% and the extra 
sources of error are all given by the N-body simulation spectra and by the estimators of the non-linear ACDM power spectrum. 

Since the fitting functions shown in the previous section are useful when applied on top of a ACDM non-linear power spectrum, 
we need to choose an estimator for the non-linear ACDM P{k). Our tests show that FrankenEmu lb], an improved version of the 
original Coyote Cosmic Emulator, performs better than the revised version of Halofit by lO, which is included in recent versions 
ofCAMB Ga. While at z = 0 both estimators work similarly well with an accuracy at the 5% level in the BAO range, Halofit 
performs much worse with increasing redshift and increasing k, as illustrated in figure At z = 0 the Cosmic Emulator shows 
a flat etTor curve for all scales up to the Nyquist frequency, below the etTor estimated for Halofit. A comparison between power 
spectrum estimators and N-body simulations has been performed also in ll73l . where they found similar results: Halofit and the 
Cosmic Emulator perform similarly in the BAO range, but the errors introduced by Halofit is above the 10% level at scales of 
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Model M7 for ;6=0.05 


Percentage Error 



Model M2 for >5=0.05 


Percentage Error 


Redshift z 





Figure 2. Error contour plot for the fitting functions of model M7 (left) and M2 (right) applied compared directly to the N-body simulations. 
We show that for the smallest coupling constant available in the simulations, j6 = 0.05, the error remains below 1% for the scales and redshifts 
we are interested in. 



Figure 3. Relative error of the fitting functions Halofit (green) and CosmicEmulator (orange) (FrankenEmu) with respect to the CoDECS 
ACDM N-body power spectra at two different redshifts (solid lines: z = 0, dashed lines: z = 1.6). While the relative error with respect to the 
CosmicEmulator remains below the 5% limit (horizontal red dashed line) for all interesting k-values, the error compared to Flalofit is already 
bigger than 10% for k > 0.5h/Mpc. Both the error of our simulation with respect to Flalofit and to the CosmicEmulator increase as a function 
of redshift. The change of trend after the Nyquist frequency (marked with the vertical grey line) can be attributed to the use of the folding 
method for the CoDECS nonlinear power spectra. 


around k « lMpc//t and z > 1. Therefore, we use the FrankenEmu as the preferred ACDM non-linear power spectrum estimator 
for our forecasting purposes. They claim to be accurate at the 1-3% percent level around the scales of interest and they have 
performed very careful resolution tests using hundreds of realizations. 

Another source of error is the sample variance error of the power spectrum when extracted from the N-body simulation. This 
depends on the number on the number n^od of independent modes available at each wavevector bin in k and its given by ll74l : 

Gp = ^J^^P{k). In figure 1^ it can be seen as a blue shaded region marked by a red and blue solid line. For large scales 
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Figure 4. Different sources of error affecting the nonlinear power spectrum. Each source of error is shown at two different redshifts, z = 0 and 
z = 1.6. The sample variance error (red solid line and blue solid line, together with their shaded regions) corresponds to the error induced by 
the limited number of available A:-modes when extracting the power spectrum; it has a sharp increase at the Nyquist frequency, since there the 
folded mesh has fewer modes available than the top level mesh, see section [Tv] for more details. The fitting error (green and orange dashed 
lines) corresponds to the intrinsic error of the fitting functions with respect to an N-body simulation using the same parameters, calculated in 
section |V] The emulator error (green and purple thick lines) is the error of the Cosmic Emulator compared to a ACDM N-body simulation 
from Codecs (these two lines correspond to the orange lines of figure]^. The relative error increases with redshift and scale and reaches 
more than 15% at the Nyquist frequency for the highest redshift. The vertical grey-dashed lines mark the scales ki^y/2 and ki^y. 


this error is considerable, but there one usually uses just linear power spectra computed by Boltzmann codes, like CAMB or 
CLASS. At k = 0.1/t/Mpc the binning error is at the percent level and decreases rapidly to negligible values, but then it increases 
again quickly at k ^ 2/i/Mpc since there the folded mesh for the power spectrum has again only few modes to sample from, as 
explained in section [TVl 

The intrinsic error of the fitting function is shown in figure]^ as a dashed green line and it remains well below the 1% level at 
all the scales of interest. 

We include all errors discussed here in our Fisher forecast analysis. It turns out that they do not affect considerably the results 
for a survey like Euclid, which will measure such a high number of galaxies, that the sampling of the clustering signal will be 
not affected by small sources of noise in the data. 


VII. FISHER MATRIX METHOD 


The Fisher matrix formalism (GSHTTl) is one of the most popular tools to forecast the outcome of an experiment, because 
of its speed and its versatility when the likelihood is approximately Gaussian (see llTSll for a method on how to improve this 
approximation). Here we apply the Fisher matrix formalism to two different probes, galaxy Clustering (GC) and weak lensing 
(WL), which are the ones that a future satellite like Euclid ||48]| will carry on board. 


A. Galaxy Clustering 

The Fisher matrix for the galaxy power spectrum has the following form 


^survey 


d\nPob,{k,ii,z) 


ddi 


dOj 


n{z)Pobs{k,IJ.,z) Y 
n{z)Pob!,{k,p,z) + l 


(9) 
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Here Vsurvey is the volume covered by the survey and contained in a redshift slice Az, n(z) is the galaxy number density as 
a function of redshift, Pohs(k^ii,z) is the observed power spectrum as a function of the redshift z, the wavenumber k and of 
/r = cos a, where a is the angle between the line of sight and the 3D-wavevector k. 

The derivatives in eq. (j^ are taken with respect to a vector of cosmological parameters, 0; = {j3^, h, log£/s, ris, COc, COh, b{z)}, 
where b{z) is the galaxy bias. When we perform the numerical derivatives of Pobs{k,iJ.,z), we have to take into account that not 
only the power spectrum and the background functions change, when the cosmological parameters are modified, but also /r and 
k are changed by a geometrical factor depending on H{z) and Da{z), due to the Alcock-Paczinsky effect ll79l . 

The observed power spectrum contains all the cosmological information about the background and the matter perturbations 
as well as corrections due to redshift-space distortions and observational uncertainties. In this work we write it in the following 
form: 


Pobs{.b^PjZ) 


Dlj{z)H{z) 

Dl{z)Hf{z) 




( 10 ) 


In the above formula, the subscript / denotes the fiducial value of each quantity, P{k,z) is the matter power spectrum, Da{z) 
is the angular diameter distance, P[{z) the Hubble function and j3d(z) is the redshift space distortion factor, which in linear 
theory is given by = f{z)/b{z), with /(z) = d\nG/d\na representing the linear growth rate of matter perturbations. The 
exponential factor represents the damping of the observed power spectrum, due to two different effects: ov. an error induced by 
spectroscopic redshift measurement errors, which translates into a damping scale Or = G^/H{z) and (Jy which is the dispersion 
of pairwise peculiar velocities which are present at non-linear scales and also introduce a damping scale in the mapping between 
real and redshift space. 


1. Including sources of error in the Fisher formalism 


In section VI we discussed several sources of error affecting the non-linear power spectrum, the intrinsic errors of the coupled 
DE fitting functions, mode-binning errors in the N-body power spectrum and the estimation and interpolation error of the ACDM 
non-linear power spectrum obtained from the CosmicEmulator. We will take these sources of error into account in our Eisher 
forecast analysis by including them as extra noise affecting the observed power spectrum. The term in square brackets in eq. 0 
corresponds to the inverse of the covariance C = P{k,z) +n{z) ^ The “noise term” n{z) ' is the number density of sampling 
points for the matter power spectrum (galaxies in a survey), which gives us an estimate of the signal-to-noise ratio we can expect 
from the forecast: for a higher number density, the power spectrum is better sampled and more information can be extracted 
from it. In order to take into account the theoretical and numerical errors on P{k,z), we decrease n(z) by a factor that contains 
the relative errors on P{k,z). In eq. 0, instead of n(z), we then have an “effective” number density: 


neff{Kz) =n{z)/{l+n(z)ap{k,z)) 


( 11 ) 


The term ap{k,z) is a scale- and redshift-dependent term which is the square root of the sum of the relative errors squared. We 
take into account all error sources which affect P{k,z) due to different reasons, as explained in section VI One of them is the 
difference between our power spectrum estimator and the N-body simulation: ap{k,z) = {Pnumericai{k,z) — Ptrue{k,z))/Ptrue{k,z). 
If dp{k,z) = 0 or it is negligible, the effective number density will be the observed one neff{z) = n{z)', otherwise, nef f{z) < n(z), 
the effective number of sampling points being reduced, together with the amount of information one can extract from the 
power spectrum. As long as n{z)P{k,z) ^ 1 for all z and k, the sampling will be always good enough to extract cosmological 
information from the power spectrum even in the presence of noise. Eor the specifications used in this work (see table [V| 
below), n{z)P{k,z) is larger than 1 in all scales of interest and therefore the theoretical and numerical error on P{k,z) does not 
have such a considerable effect on the parameter estimation, as one would expect naively. We test the inclusion of the effective 
number density ngff{z) on the Eisher forecast analysis and we find that the relative l-cr marginalised errors on each cosmological 
parameter are between 8% and 15% higher when using the estimated uncertainties on the power spectrum. 


2. Systematic bias on cosmological parameters 

In this section we will quantify the effect of the systematic errors due to the uncertainties on the non-linear power spectrum. 
We will show how big this systematic bias would be, if we used for our forecasts a power spectrum which is not the “correct” 
one, for example one obtained from a single N-body realization. Eor the following discussion on systematic biases, we will use 
the expressions derived in the Appendix B of lISOl . 





10 


The linear bias on a cosmological parameter 5 9i due to the bias 5 xj/i in a parameter of the model which we assume fixed and 
cannot be measured is given by: 


59i = - 



( 12 ) 


In our case we will have only one systematic parameter xf/, which controls the difference between the “true” power spectrum 
Pfrue and and our simulated power spectrum P„um- 


Py — ^Pfuim “f (1 ^}Ptrue 


\j/ can vary continously so that for y/ = 1 we recover P„um, while for y/ = 0 we obtain Pmie- We can define the relative difference 
between P,n,e and as: 


<yp{k,z) — 


^num {k,z) Ptrue {k,z) 


Ptrue{k^z) 


The in eqnll^above is simply the usual Fisher matrix: 


F®® = itr 


(c-^dl^cc-^dfc 


(13) 


while the pseudo-Fisher matrix between measured and assumed parameters F®*^ is: 


Fij'f' = itr (c-^dfCC-^dJC^ (14) 

which for one systematic parameter only, is just a column vector. 

In the case of galaxy Clustering we will compute F®*^ in the following way, using the fact that for y/ = 1, C — Pnum{k,z) + 
n-^{z) and Fy, 1^=1 = P„um- 




dkk^ 


^effiz)Pnumik<^z) 
XXeff{z)Pni(m{k^Z^ 1 


i{k,z) 


dP^ 

dyf 


dP„, 


V'=1 


d9i 


(15) 


in this step we have assumed that we have no systematic parameters affecting the galaxy number density n{z) and therefore, 
its derivative w.r.t y/ vanishes. Also, for notational simplicity we left out the integral over /t and the form of the observed power 
spectrum. 

We have then dy, 


1 true + Pmim — Ptrue(yp{k,z) and we just need to perform eqn. 
systematic biases on the cosmological parameters. 


■'\j/P\j/ — Ptr 


15 


with eqn. 
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in order to obtain the 


In the following table IV we present the results on the systematic bias, for a standard ACDM forecast, for different maximum 
k coverages, up to a maximum k of l.lh/Mpc. We will regard as a “true” power spectrum Ptrue, the one obtained by the Cosmic 
(Franken) Emulator ||6l, since they have performed a careful analysis of resolution effects using a large set of simulations and 
claim to be accurate for k < l.Oh/Mpc at the 1% percent level compared to state of the art N-body simulations. On the other 
hand, the numerical “biased” power spectrum F„„„„ is the one obtained from the CoDECS ACDM run, which consists on only 
one realization. We have left out the CDE coupling parameter j3, since in that case we do not have any other prediction in the 
non-linear regime to compare with. 


Parameter 

h 

Injs^ 

ns 

COb 

(Oc 

fiducial 

0.7036 

2.42 

0.966 

0.04503 

0.2256 

kmax = k]Vy/2 = l.lh/Mpc 

syst. bias with 8yf = 1 

-0.0016 

-0.15 

0.061 

0.0028 

-0.0031 

kmax = 0.35h/Mpc 






syst. bias with 8yf = 1 

-0.0094 

-0.11 

0.045 

0.0021 

-0.0039 

kmax = 0.15h/Mpc 






syst. bias with 8yf = 1 

-0.0032 

-0.04 

0.018 

0.00026 

-0.0024 


Table IV. Systematic bias on the ACDM cosmological parameters evaluated at our fiducial model. The systematic bias is higher when probing 
smaller scales, where uncertainties in the non-linear power spectrum are higher. This highlights the importance of modelling correctly the 
non-linear power spectrum in order to analyze data. Using the wrong non-linear power spectrum produces statistical errors which are larger or 
of the same order as the statistical results. 
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The results of table IV show how important it is to model accurately the non-linear power spectrum in order to make forecasts 
and to analyze the upcoming data of large scale structure surveys like Euclid ll48l . The systematic errors on the cosmological 
parameters can be bigger than the statistical errors (compare to table [V^ and figure |^in the results section below). This is the 
case if, as in our example scenario, we would use a non-linear power spectrum that is inaccurate by about 10-15% in the non¬ 
linear regime at higher redshifts (which was shown in figure]^. Therefore, it is well justified to choose for our Fisher forecasts 
the Cosmic Emulator as the “true” ACDM non-linear power spectrum estimator, since this is the most accurate predictor up to 
date. It would still be interesting to know how robust is the signal of the coupling parameter j3 with respect to changes in the 
other parameters or in the ACDM non-linear prediction, but as long as we do not have better and faster semianalytic methods 
applicable to general models of dark energy, the estimation of systematic biases of extra parameters is an impossible task. 


3. Choice of the k„^ax integration limit in the Fisher formalism 

There are at least three ways of setting the maximum k mode used in the Fisher forecast integration (eqn. |^. One common 
choice is to set a hard cut in k at all redshifts, and depending if one wants to include or not non-linear effects, this cut can be 
taken at linear scales, smaller than k = 0.1/z/Mpc, or at non-linear scales k > 0.1/z/Mpc. In the latter case if one is using a 
power spectrum calculated in linear theory, one needs to to use some Lagrangian damping correction, as introduced originally in 
Gguni, in order to take into account broadening effects on the BAO peak induced by the non-linear evolution of densities and 
velocities. Another option to cut off the power spectrum is to demand that the variance (y%{z',k) stays below a specified value at 
each redshift, therefore implicitly changing the cutoff in k as a function of z. An usual choice for this is 0%{z',kcut) = 0.25, as 
was done previously in Il82l . Since we are assuming to have a knowledge of the non-linear power spectrum up to the Nyquist 
frequency, we implement for our forecasts a hard-cut method, at kf^y and at %y/2, without the need of any Lagragian damping 
correction. However, to be conservative, we cite as our main results the ones in which the cut is performed at kf^y/2 in order 
to eliminate any possible unknown contribution from the numerical high-frequency noise entering the estimation of P{k) at the 
Nyquist frequency (see e.g. Il^l83]l for similar prescriptions on where to cut the non-linear power spectrum). 


B. Weak Lensing 

For the weak lensing power spectrum the Fisher matrix is just a sum over all possible correlations at different redshift bins 
ITT?! . namely : 


^ap — Jsky 2^ 2^ ^ ^ .-t ^ ^yr 

£ ij,k,m 


2 dda P d da 


(16) 


With the power spectrum of the shear field given by: 


9 Wiiz)Wjiz)HAz)Dliz) 




PmA/r{z)) 


(17) 


where P,n is the matter power spectrum discussed above, the indices i, j stand for each of the redshift bins and the window 
functions are given by: 


W{z) = l^ dz(^l-^)n(f) 


(18) 


where the normalized galaxy distribution function is: 


n{z) = z^exp (-{z/zoA^^^ (19) 

Here the median redshift and zo are related by Zmed ~ 1.412zo. The corresponding covariance matrix has the following form: 


Qji£)=Ptji£) + 5ij^„,nA+Kijie) 


( 20 ) 


where jint is the intrinsic galaxy ellipticity and the shot noise term contains 
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rii = 3600 



2 

uqI-jY 


( 21 ) 


with hq the total number of galaxies per arcmin^ assuming that the redshift bins have been chosen such that each contain the 
same amount of galaxies (equi-populated bins). The matrix Kij{£) is a diagonal “cutoff” matrix whose entries increase to very 
high values at the scales where the power spectrum P{k) has to be cut to avoid the inclusion of numerical noise errors. Since 
Kij{£) depends on the multipole £, which itself depends on k and z through: £ = kr{z), each of the entries Ka increase with the 
multipole in a different way. For our purposes we chose a very sharp cutoff, such that multipoles containing wavenumbers k 
beyond the Nyquist frequency (or half of the Nyquist frequency, for our reference case) at the center of the redshift bin i, do not 
contribute to the Fisher matrix at all after the corresponding £cut (specified in table VIII below). 


C. Adding non-linear corrections to the power spectrum 


We are interested in a Fisher forecast that includes information from non-linearities. In this case we cannot separate the 
power spectrum, into a power spectrum at redshift zero multiplied by the square of the normalized growth factor, but we need 
to evaluate directly the non-linear power spectra P{k,z) at different redshifts. A full non-linear correction for the redshift space 
distortions would be also desirable, but since the modelling of that effect is yet to be understood in general cases, we use as a 
first approximation an exponential damping of the form exp{—k^jj.^ay), where Gv is the pairwise peculiar velocity of galaxies 
induced by non-linearities in the matter and velocity power spectrum. This is the hrst term of a set of corrections that can 
be applied to the Kaiser formula ll84l for the clustering in redshift space (see e.g. ia, IMl, El,®). We use a value of 
(Jy = 300km/s which is an approximate and conservative value based on the estimations by lf85l . in which the authors test a 
variety of redshift-space distortion models. 

As already mentioned, the damping term, introduced by G6l, which should correct the linear P{k), especially the position of 
the BAO peaks, for non-linearities, is not included here, since we assume that we have a complete model of the non-linear power 
spectrum and therefore all possible corrections are already included in our fitting functions and power spectrum emulators. 

In order to implement our model of coupled Dark Energy, we use for H{z) and D{z) tables precomputed using a modihed 
version of CAMB, that includes the exponential and inverse-power law potentials for coupled dark energy llJTlIS^ . We calculate 
these tables for each of the parameters 0, ± £, where e = 0.03. We do the same for the linear perturbation quantity G{z), whose 
logarithmic derivative with respect to ln(a), known as the growth rate /(z) enters the redshift space distortion term in 10 The 
background quantities H{z) and D{z) are important in the Fisher forecast for the Alcock-Paczynski geometrical term, which we 
take into account in the observed power spectrum. 

The full non-linear power spectrum we use in our method is then obtained as follows. 


• At linear scales k < 0.1/i/Mpc the linear power spectrum is obtained from our modified version of CAMB which includes 
the effect of the DM-DE coupling. 

• At non-linear scales, k > 0.1/i/Mpc, we use a combination of the power spectrum calculated with the cosmic emulator 
FmnkenEmu and our fitting formulas for CDE that account for the non-linear dynamics in presence of a fifth force. 

• The matching at k « O.lh/Mpc is performed using different interpolation methods, either averaging on both sides and 
smoothing out or allowing for a small discontinuity. The matching point is left out of the Eisher integral and we check for 
this effect, showing that the different methods have a negligible effect (less than 2 percent) on the hnal absolute errors on 
the parameters. 


In order to be conservative in terms of numerical etTors and noise, we cut the power spectra at half of the Nyquist frequency 
(we also test and compare with a cutoff at kffy) and we include all sources of errors specified in section VIIA 1 , as effective 
noise terms into our Eisher estimation. 


VIII. RESULTS 


We now present results for galaxy clustering and weak lensing, using the specifications for a Euclid-like survey as described 
in tables [V| and |VII| below. We use the Eisher formalism described in section VII to forecast the errors in the cosmological 
parameters, using information from the non-linear power spectrum for coupled Dark Energy models, as obtained with the 
procedure described in section|IV]together with the fitting functions obtained in section|V] To make our estimation more realistic, 
we also take into account all sources of error and systematics for the power spectrum which were discussed in section VI As a 
way of testing our improvement on the parameter estimation, we also perform two extra Eisher forecasts, first using only linear 
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Parameter 

Value 

Description 

^survey 

15000 deg2 

Survey area in the sky 


0.001 

Photometric redshift error 


300 km/s 

Fiducial pairwise velocity dispersion 

Az 

0.2 

Redshift bin width 


{0.6, 2.0} 

Min. and max. limits for redshift bins 


Table V. Specifications for the Fisher Matrix of an Euclid-like galaxy clustering survey. 


Parameter 


h 

lO'^ss/s 

ns 

Q.h 

klc 

fiducial 

0.0025 

0J036 

2.42 

0.966 

0.04503 

0.2256 

l^max — 







abs. error 

0.000346 

0.00160 

0.01855 

0.00267 

0.00084 

0.00088 

relative error 

14% 

0.23% 

0.77% 

0.28% 

1.9% 

0.39% 

kmax = 







abs. error 

0.000305 

0.00151 

0.01808 

0.00249 

0.00081 

0.00085 

relative error 

12% 

0.22% 

0.75% 

0.26% 

1.8% 

0.38% 


Table VI. 1-CT fully marginalized errors on the cosmological parameters of the model for a galaxy clustering Fisher forecast. Two cases are 
presented, a hard cutoff of the power spectrum at kmax = and a hard cut at ky^y. 


power spectra for the CDE model and then correcting these linear P{k) with the latest version of Halofit |[3l, which was designed 
for ACDM only. 

The fiducial parameters are cOc = 0.1117, (Ob = 0.0223, = 0.966, = —19.8395, h = 0.7036, = 0.0025, which 

are consistent with WMAP7 results (see table [n|. The fiducial values for the galaxy bias b{z) used for the GC probe are 
taken from the Euclid specifications (see El |83) For our final results, we convert these parameters into the set pi = 
h, ns, Q.C, ^b} marginalizing over the bias b{z) (for the GC case) and using a Jacobian transformation to convert 

into the new set of parameters, which is allowed by the Eisher matrix formalism. We choose to forecast the error on the square 
of the coupling parameter j3^, because this is the quantity entering the modified gravitational Newton constant Geff in the limit 
of linear perturbations (cf eq. [7]l, therefore giving the strength of the “fifth force”. The corresponding fiducial value for the 
coupling, P = 0.05, is still compatible with recent limits set by analyzing the data from the Planck Satellite (see ElIlQl). 


A. Galaxy clustering 


Table [V] shows the specifications of a Euclid-like survey, which are used in our Eisher forecast. While Euclid specifications 
use 14 redshift bins of a width Az = 0.1 (see ll48l . table 3 in that work), we use only 6 bins of a width of 0.2. We check in the 
case of ACDM, that this re-binning ( done using the specified number of galaxies and the corresponding comoving volumes in 
our cosm ology ) has a very small effect (of a few percent) on the estimation of the l-cr errors on the parameters. 


In table 


VI 


’ we show the fully marginalized 1-ff errors on the final cosmological parameters pj = {p^, h, 10^, 2 ^^, ris, Q.c, 
performing the non-linear power spectrum cut-off at two different scales: k^ax = kNy/2 and k,nax = kiyy. As explained above, 
we take as a reference result the one corresponding to a cutoff at ky/y/l. The gain in constraining the parameter when going 
from kf^yjl to ky^y is of two percent points in the relative errors, while for the other cosmological parameters the improvement is 
negligible. 


1. Variation of the kmax integration limit 


We now test the gain in information obtained by including progressively more non-linear wavemodes k into the Eisher inte¬ 
gration. We perform the same Eisher forecast, each time increasing the maximum mode k^ax at which the integration is cut off. 
In figure|^we show how the l-CJ fully marginalized error on the cosmological parameters p, changes with an increase of kmax- 
The error decreases steadily with an increase of kmax, where the biggest gain is achieved when going from linear {k « O.l/z/Mpc) 
to mildly non-linear (k « 0.3/i/Mpc) scales. Eor parameters like h and Q.b, an approximate plateau is reached already before 
kmax ~ 10, while for P^ there is still a considerable gain in parameter constraints when going into smaller scales, even beyond 
kyiy/2 (consistent with table VI above). This happens, qualitatively, because at small scales we have a well-defined characteristic 
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Figure 5. Change of the 1-CT fully marginalized error on the set of cosmological parameters pi for a galaxy clustering Fisher forecast, as a 
function of the maximum mode k^ax used as a cutoff in the Fisher matrix integration. When increasing the maximum k,„ax the errors on the 
parameters get steadily smaller, especially when going from linear {k ~ O.l/r/Mpc) to mildly non-linear {k pa 0.3/i/Mpc) scales. The vertical 
dashed grey lines, mark the half and the full Nyquist frequencies. 
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Figure 6. Marginalized confidence contour regions 1,2,3-0' for jS^ and h. The plots correspond to a comparison between linear (green-dashed 
lines) and non-linear (red-solid lines) scales for GC (left) and WL (right) using the scale cutoff at kj^y/l. 


signal coming from non-linear interactions including the fifth force 'which itself involves the p coupling, while the information 
on the backgorund cosmological parameters gets washed out (c.f. section V C above). 
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Parameter 

Value 

Description 

I sky 

0.364 

Observed sky fraction 

Tint 

0.22 

Intrinsic alignment 

ne 

30 

Galaxy density per arcmin^ 

^pz 

0.05 

Photometric redshift error 


Table VII. Specifications for the Fisher Matrix of an Euclid-like weak lensing survey. 


Parameter 

Value 

Description 

.A' 

^peak 

^cut,k] 

^CUt,k2 

Zrange 

6 

{0.59,0.75,0.90,1.06,1.28,1.57} 
{0.22,0.23,0.25,0.27,0.32,0.51} 
{1686,2070,2410,2753,3155,4344} 
{3372,4141,4820,5506,6311,8689} 
0.5<z<3 

Number of redshift bins 

z-position of peak of the bin 
full width at half maximum of the peak 
cutoff in multipole 1 at the center of each bin for k,„ax = kjgy/2 
cutoff in multipole 1 at the center of each bin for k,„ax = k^y 
Total range in redshift of each bin 


Table VIII. Redshift bins specifications for an Euclid-like weak lensing survey using equipopulated redshift bins in the range 0 < z < 3 and 
the corresponding values for the cutoff applied in the multipoles (. at each redshift bin, for two different cases: scales larger than k\ = kj^y/l 
and scales larger than ^2 = 


B. Weak lensing 


Table VII show the specifications for a weak lensing probe in an Euclid-like survey. The redshift bins are chosen in such a 
way that they contain an equal number of galaxies (equipopulated bins). The bins are then given by: 


ni{z) = :^n{z) 


Erf 


/ Zi+l -z 


^ -Erf 


Zi — Z 




( 22 ) 


where f,- are the values of the bin intervals in the range Zrange = 0-5 < z < 3 chosen such that for each interval the integral over 
the galaxy distribution function n{z) (eqn. 19 1 is equal. The resulting peaks of the bins and their full width at half maximum are 
specified in table [VIlI 


Analogously to the galaxy clustering case, we show in table|I^the fully marginalised 1-ff errors on the parameters p,-. The 
sum over multipoles f in eq. ((T6]l is performed from = 5 up to a maximum of = 20000, but as explained in section 
VIIB we perform a cutoff at each redshift bin, so that no scales in the non-linear power spectrum beyond the half of the Nyquist 
(for our reference case) or beyond the Nyquist frequency (for our second case) contribute to the WL signal. The values of these 
cutoffs, (.cut for the two different cases k\ = kNy/2 and k 2 = k^y are listed in table VIII In contrast to the GC case, going from 
kfjyjl to kf^y in a WL survey does bring a noticeable improvement on the estimation of parameters. 

The dependence of the l-CJ fully marginalized error on each parameter p, with respect to (max is shown in figure]^ When 
using just linear power spectra information, the error on the parameters does not improve if one increases the scale (max, while 


Parameter 

J3' 

h 

10 Vi 

ns 

Q.b 

O.C 

fiducial 

0.0025 

0.7036 

2.42 

0.966 

0.04503 

0.2256 

4ut.kN,/2 







abs. error 

0.000125 

0.00835 

0.112 

0.0105 

0.0032 

0.0046 

relative errror 

5.0% 

1.2% 

4.6% 

1.1% 

7.1% 

2.0% 

^CUt.kNv 







abs. error 

0.000097 

0.0068 

0.058 

0.0085 

0.0022 

0.0032 

relative error 

3.9% 

0.97% 

2.4% 

0.88% 

4.9% 

1.4% 


Table IX. I-c fully marginalized errors on the cosmological parameters of the model for a weak lensing Fisher forecast. Two cases are 
presented, a redshift-dependent cutoff kyj! ® cutoff icut^ky^ corresponding to cutting off the non-linear power spectrum at the half and 
at the full Nyquist frequency (analogously to the GC case) as explained in the main text. As opposed to the GC case, going into smaller scales 
in a WL survey does bring a noticeable improvement. 
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Figure 7. Variation in the l-CT fully marginalized errors for each cosmological parameter pi as a function of the maximum multipole used in 
the weak lensing Fisher forecast. Increasing the observed multipole range improves considerably the determination of a parameter, especially 
on the coupling j3^ (red solid line) and the initial amplitude of perturbations 10®(dashed green line). The double lines corresponding to each 
parameter, show how the error changes if a cut in the matter power spechum P{k) is performed at kfjy (lower line) instead of a cut at kj^y/l 
(upper line, respectively). The vertical dashed grey lines mark the tcut at the peak of the first redshift bins for the cases ksy/'i. ^cut.kny- 


in the case where non-linear information is used, increasing the maximum multipole l,nax improves considerably the l-(7 error 
on the parameters, especially on the coupling and on the initial amplitude of scalar fluctuations 10®As. This is due to the 
fact that the extra signal on the coupling coming from the non-linear part of the power spectrum, the so called “bump”, greatly 
enhances the constraints on the parameter estimation. The double lines corresponding to each parameter in figure]^ show how 
the error estimation is changed if scales up to kf^y are included (lower line) compared to the upper line where only scales up to 
kf^y/2 contribute. At small i both lines are on top of each other and only start diverging at around £ = 2000, when the extra 
amount of information contained in highly non-linear scales starts becoming important. The most significant gains occur again 
on the parameters j5^ and lO^A^. 


C. Combined results 

In figure]^ we show the 1-, 2- and 3-(7 confidence contours from the Fisher forecast for WL and GC. These confidence 
regions for each pair of parameters are obtained after marginalizing over all the other parameters. As it can be seen, some 
degeneracies are broken when combining the confidence ellipses from two different observations, for example in the case of 
the plane Q.ti, j3^. Other parameter combinations, as show the same orientation of the ellipses for WL and GC, so that 

the combination of both probes does not help to disentangle the degeneracies. While GC constraints much better the usual 
parameters {h, 10®Mi, Oc O*}, WL constrains the coupling parameter much better which can be seen in the vertical 
orientation of the ellipses that correspond to )3^. Therefore combining the observations on GC and WL, as a future survey like 
Euclid will do, is a powerful way of constraining degenerate parameters in cosmology. 

In tablewe cite the l-CJ fully-marginalised errors on the parameters pj for three different cases: a) using only linear CDE 
power spectra computed from our modified version of CAMB; b) applying a non-linear correction to these linear CDE power 
spectra using the latest version of Haloflt from 13, which was designed for ACDM-only; c) using the full coupled DE non-linear 
power spectra computed with our fitting functions following the procedure explained in section [VIIC| 

This shows the value of using the N-body-calibrated fitting functions on the coupling Using the proper j3-dependent non¬ 
linear correction instead of the standard Haloflt correction, the constraints on improve by more than an order of magnitude for 
WL and by a factor of order three for GC. This makes very clear the importance of applying non-linear corrections that depend 
on the parameter to be tested. 
























17 




Figure 8. Marginalized confidence contour regions (1,2,3-0') for all cosmological parameters considered in this model. The blue dashed 
lines correspond to the WL Fisher forecast, while the green dashed lines correspond to the GC Fisher forecast both in our reference case. The 
red solid lines correspond to the combined Fisher matrix forecast. For combinations of the parameter /3^, WL and GC have similar figures 
of merit, but different orientations, while for other combinations of cosmological parameters, the estimation is dominated by the GC Fisher 
matrix estimation. 
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Parameter 

|S" 

h 

lOVs 


fife 

flc 

fiducial 

0.0025 

0.7036 

2.42 

0.966 

0.04503 

0.2256 

WL: 1-cT abs. error, using: 







linear CDE 

0.0189 

0.040 

0.221 

0.0139 

0.0062 

0.0127 

linear CDE-l-Halofit 

0.0184 

0.044 

0.256 

0.0109 

0.0066 

0.0079 

non-linear CDE fitting functions 

0.000125 

0.00835 

0.112 

0.0105 

0.0032 

0.0046 

GC: 1-<T abs. error, using: 







linear CDE 

0.0038 

0.0097 

0.117 

0.0176 

0.0021 

0.0055 

linear CDE-l-Halofit 

0.0011 

0.0029 

0.024 

0.0023 

0.0007 

0.0006 

non-linear CDE fitting functions 

comb. WL-i-GC: l-cr error, using: 

0.00035 

0.0016 

0.018 

0.0027 

0.0008 

0.0009 

non-linear CDE fitting functions (abs.) 

0.000084 

0.0010 

0.0169 

0.00251 

0.00072 

0.00080 

non-linear CDE fitting functions (rel.) 

3.4% 

0.16% 

0.7% 

0.26% 

1.6% 

0.35% 


Table X. l-CT fully marginalized errors on the cosmological parameters for WL, GC and the combined Fisher matrix WL+GC, using three 
different power spectra. Linear CDE: Using only information from the linear power spectrum for the CDE model up to a scale of k = 
0.1/i/Mpc. Linear CDE+Haloflt: Using the linear power spectrum for CDE plus a non-linear correction using the latest Flalofit from a 
Non-linear CDE fitting functions: Using the fully non-linear power spectra for CDE obtained from the fitting functions and the emulator 
as explained in|V] which we regard as the most reliable description of the model in this range of scales. In all these cases we are using our 
reference Fisher forecasts corresponding to the cutoff at kfjy/1. 


IX. CONCLUSIONS 


The goal of this paper is to exploit the cosmological information contained in the non-linear regime in order to improve 
parameter estimation from future large-scale observations. The main obstacle along this road is that we have accurate non-linear 
corrections for the matter power spectrum only for ACDM and a few other relatively simple variants, but not for the large variety 
of modifed gravity models that have been proposed in recent years. 

The first part of this paper has then been devoted to the task of finding corrections to the linear power spectrum in the range 
of k « 0.1 — Ih/Mpc for a selected class of modified gravity models, namely coupled dark energy. This model is indeed one 
of the simplest possible extensions of Einstein’s gravity and depends entirely on a single parameter, the coupling constant j3 (in 
addition to the standard ones). Employing the CoDECS suite of simulations we build different fitting function models, such 
that when multiplied by the ACDM non-linear power spectrum (we use the estimator provided in ref. ) they reproduce the 
N-body results to an accuracy of 1%, for scales k between 0.1 and 5 h/Mpc and a range in z, between 0 and 1.8. To achieve 
this accuracy in the fitting functions we need to perform a careful extraction and interpolation of the power spectrum from the 
simulation mesh. 

The accurate fitting functions have been then employed to extend the regime of validity of the forecasts for future experiments. 
We focused on a Euclid-like survey that includes weak lensing and redshift-space distortions (galaxy clustering) and predicted 
the constraints in the cosmological parameters, with particular emphasis on the dark matter-dark energy coupling j3. We find 
that j3 is better constrained by weak lensing than by galaxy clustering (contrary to all the other standard parameters). We find 
that the extension into non-linear scales improves the constraints by more than an order of magnitude compared to previous 
results using only linear power spectra, but also by more than an order of magnitude in WL and a factor of three in GC compared 
to using a wrong ACDM Halofit non-linear correction. We also show that using the wrong non-linear power spectrum, can 
bias systematically the estimation of errors on the cosmological parameters, yielding systematic errors of the same order of 
magnitude as the statistical ones. This makes very clear that it is important to include the proper non-linear corrections to the 
parameter to be tested, especially for models beyond ACDM in which the small-scale gravitational dynamics are modified. 

To make our forecast more realistic, we take into account all known sources of error entering the estimation of the power 
spectra and the fitting functions in the way of a reduced effective number density of galaxies and then perform a conservative cut 
of the power spectrum at half of the simulation Nyquist frequency, to avoid other sources of unknown numerical noise affecting 
the results. In the case of GC we include also a first approximation to the correction to redshift space distortions, caused by 
peculiar pairwise velocities at non-linear scales. We find that a space probe like Euclid will be able to constrain the coupling 
parameter around the fiducial value 0.0025 at l-(7 with a relative accuracy of 14% when using weak lensing alone, 5% when 
using only galaxy clustering and at 3.4% when combining both probes. 

It is interesting to note that the most stringest constraint we obtain amounts to Aj3^ « 8 • 10^^; this level of precision on the 
dark matter-dark energy coupling is not far from the current best limits reached with Solar System observations on a coupling to 
baryons ll52]| . which can be translated in our notation as < 2 • 10^^ at l-cr. 
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Appendix A: Parameters of the fitting models M2 and M7 
Model M2 and M7 have each 5 coefficients, which are third order polynomial function of j3 and z: 

at (j3, z) = qa + qap + qa z + q^pz 

+ qi5P^+qi6Z^ + qn z^+qisli^ (A1) 

Since we have 5 coefficients and 8 free parameters for each coefficient, we have a total of 40 free parameters. They 
can be represented in a matrix qtj by: A, = QijBj, where the vector A, = ai(l5,z) is the vector of coefficients and Bj — 

(1, j3, z, j3z, z^, z^) the vector of independent variables and the indices are defined in the ranges: i = 1, j = 1, 

The corresponding matrix qtj for model M7 is: 

/0.00689754 -0.200001 -0.0112943 0.157222 -1.40171 0.0129681 -3.99345 -0.00629434\ 
-0.0706221 1.45187 0.236402 -1.26467 9.13266 -0.300832 9.81754 0.120649 

qij{Ml)= 0.069102 -1.44396 -0.24268 1.26135 -9.23052 0.311995 -12.1568 -0.12561 (A2) 

1.19797 0.632936 -1.74767 0.920634 -8.47056 1.46434 42.7015 -0.471948 

\ 0.684838 1.15722 -2.15405 3.93688 -5.29011 3.42683 -43.9295 -1.18691 

and the coefficient matrix for model M2 is: 

/0.00301472 -0.0961879 -0.00897697 0.0389941 -0.860642 0.0131329 -4.73786 -0.00552858\ 

-0.028584 0.615003 0.111671 -0.371008 4.94379 -0.138279 20.3258 0.0512362 

qij{M2)= 0.0132728 -0.345125 -0.0701132 0.115389 -3.11777 0.102871 -16.9949 -0.0408223 (A3) 

2.91757 1.03545 -5.16595 1.69976 -22.1649 4.33453 120.513 -1.35074 

\ 0.867566 2.92363 -1.14174 3.4356 -2.96193 1.91413 -113.348 -0.684269 

Applying a singular value decomposition on both matrices, shows that they have a very similar structure, therefore since the 
sigmoidal functions M2 and M7 have a related functional form, the matrices describe roughly the same information. 



